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^solved  .in  either  the  forward  or  the  adjoint  mode.< 


Time-dependent  fluxes,  and  flux  functionals  such  as  dose, 
heating,  count  rates,  etc.,  are  calculated  as.  function  of 
energy,  time  and  position.  Multiple  scoring  regions  are 
permitted  and  these  may  be  either  finite  volume  regions  or 
point  detectors  or  both.  Scores  at  point  detectors  are 
made  by  a special  "bounded  estimation”  procedure  which 
elimihates  unbounded  scores  due  to  the  usual  inverse-squar 
estimation' process.  Other  scores  of  interest,  e.g.,  # 
collision  and  absorption  densities,  etc.,  are  also  made. 

A special  feature  of  SAM-CE  is  its  use  of  the  '"combina- 
torial geometry"^  technique  which  affords  the  user  geometrii 
capabilities  exceeding  those  available  with  other  commonly 
used  geometric  packages.  An  automatic  geometry  checker  is 
included.  . 

All  nuclear  interaction  cross  section  data  (derived  from  • 
the " ENDF  IV  libraries)  are  tabulated  in  point  energy 
meshes.  The  energy  meshes  for  neutrons  are  internally 
derived,  based  on  built-in  convergence  criteria  and  user- 
supplied  tolerances.  Tabulated  neutron  data  for  each 
distinct  nuclide  are  in  unique  and  appropriate  energy 
meshes.  Both  resolved  and  unresolved  resonance  parameters 
from  ENDF  data  files  are  treated  automatically'  and  ex- 
tremely precise  and  detailed  descriptions  of  cross  section 
behavior  is  permitted.  Such  treatment  avoids  the  ambi- 
guities usually  associated  with  multi-group  codes,  which 
use  flux-averaged  cross  sections  based  on  assumed  flux 
distributions  which  may  not  be  appropriate. 

By  use  of  the  "band"  feature  of  the  code,  which  automati- 
cally splits  cross  section  <iata  into  two  or  more  energy 
ranges  to  be  treated  one  at  a time,  SAM-CE  affords  one 
the  ability  to  consider  many  nuclides,  in  a given  configur- 
ation, each  being- described  in  much  detail. 

SAM-CE  also  provides  the  user  with  the  opportunity  to 
employ  energy,  region  and  angular  importance  sampling. 

An  extensive  library  of  processed  cross  sections  for  use  ir 
SAM-CE  is  available  at  the  Radiation  Shieding  Information 
Center  (RSXC) , ORNL,  Oak  Ridge,  Tennessee. 
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SECTION  2 - PROGRAM  SAM-X 


2.1  General  Description  of  the  SAM-X  Program 

With  the  establishment  of  the  Cross  Section  Evaluation  Center 
at  Brookhaven  National  Laboratory,  detailed  and  up-to-date  evalu- 
ated cross  section  information,  in  the  form  of  the  ENDF  files,  has 
become  available.  In  order  to  make  use  of  this  wealth  of  accurate 
information  in  the  SAM-CE  System,  the  SAM-X  program  which  generates 
processed  cross  section  data  tapes  for  later  use  in  the  SAM-F  and 
SAM-A  Monte  Carlo  codes,  was  written. 

SAM-X  is  designed  to  process  ENDF  neutron  and  gamma  ray  cross 
section  and  gamma  ray  production  data  files.  The  output  of  SAM-X 
is  a neutron  element  data  tape  (NEDT) * , a gamma  ray  production  data 
tape  (GPDT) * , and  a gamma  ray  element  data  tape  (GEDT)*,  which  are 
subsequently  used  as  input  to  the  SAM-F  and  SAM-A  Monte  Carlo  trans- 
port codes. 

SAM-X  is  a program  with  an  overlay  structure.  The  code  com- 
prises a small  driver  (main  overlay),  and  five  processors  (primary 
overlays):  NUTRON,  PEND,  WEED,  GAMMA  and  BCDEAN.**  Program  NUTRON 

processes  the  ENDF  neutron  files.  Program  PEND  (Production  from 
Evaluated  Nuclear  Data)  processes  the  ENDF  photon  production  files. 
Program  WEED  (Weed-out  Extra  Energy-dependent  Data) , which  comple- 
ments PEND,  is  a data  reduction  code.  Program  GAMMA  processes  the 
ENDF  gamma  ray  cross  section  files.  Program  BCDEAN  converts  the 
normal  binary  mode  output  of  SAM-X  to  BCD  mode  for  inter-facility 


* In  the  descriptions  that  follow,  EDT  (element  data  tape)  will 
generally  refer  to  either  an  NEDT,  a GPDT,  or  a GEDT,  where  the 
specific  reference  intended  is  clear  from  context. 

**  Referred  to  as  NUTRON,  PEND,  WEED,  GAMMA  and  BCDEAN  throughout 
Section  2,  these  primary  overlay  programs  are  actually  coaed  as 
PROGRAM  OVERIO,  OVER20,  OVER30,  OVER40  and  OVER50,  respectively, 
as  convenient  flags  in  switching  from  CDC  to  IBM  versions. 

15 


data  reduction.  But  WEED  may  be  utilized  to  reprocess  a GPDT 


which  has  already  been  "WEED-ed",  if  so  required.  In  addition, 
the  internally  set  default  criteria  may  be  superseded  by  user- 
specified  criteria. 

Program  GAMMA  reads  the  ENDF  gamma  ray  cross  section  file 


(File  23).  Since  both  Monte  Carlo  codes,  SAM-F  and  SAM-A,  do  not 


consider  coherent  scattering  at  the  present  time,  the  total  cross 
sections  (MT«501)  given  by  ENDF  are  modified  by  subtracting  out 
the  coherent  scattering  (MT=502) . GAMMA  then  treats  incoherent 
scattering  (MT-504)  and  re-tabulates  it  in  the  same  energy  mesh  as 
the  modified  total  cross  sections.  Note  that  the  Monte  Carlo  pro- 
grams subsequently  use  the  classic  Klein-Nishina  distribution  for 
scattering  so  that  an  inconsistency  develops  in  the  low  energy 
(x-ray)  range.  The  organization  of  a GEDT  for  an  element  is  des- 
cribed in  Appendix  B„ 

The  machine  dependence  of  the  binary  EDT  produced  by  the  first 
four  primary  overlays  motivated  the  development  of  the  post  pro- 
cessing program,  BCDEAN.  This  fifth  primary  overlay  serves  as  a 
two  way  EDT  converter,  i.e.,  binary-to-BCD  and  vice  versa.  The 
structure  of  a typical  BCD  EDT  is  described  in  Appendix  N. 

The  user  may  now  proceed  to  Section  2.5,  2.6  and  2.7  which 
contain  descriptions  of  the  input,  additional  notes  to  the  user. 


and  tape  and  file  utilization,  respectively.  However,  those  de- 
siring a more  thorough  background  of  SAM-X  may  wish  to  read  Section 
2.2,  2.3  and  2.4  which  contain  a general  description  of  the  ENDF 
System,  the  theory  behind  the  SAM-X  code,  and  the  program  informa- 
tion flow,  respectively. 


m ( 


2.5  SAM-X  Processed  Cross  Section  Data 

As  an  aid  to  the  user,  files  of  SAM-X  processed  ENDF  data  are 

available  through  the  Radiation  Shielding  Information  Center  (RSIC) 

at  Oak  Ridge  National  Laboratory.  Table  2.6  shows  the  nuclides  for 

which  SAM-X  neutron  data  files  are  available.  The  right  hand  column 

of  Table  2.6  shows  for  which  of  these  nuclides  gamma  ray  production 

* 

data  files  are  also  available. 

Note  that  Table  2.6  is  as  of  31  October  1975.  It  is  antici- 
pated that  the  number  of  processed  files  will  be  increased  in  the 
future. 


A complete  file,  for  all  elements  (Z»l  to  94)  is  available  for 
gamma  ray  cross  sections. 
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1002 
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Mod 

0 

Li-6 
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1271 
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3007 
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* 

Be-9 
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Pb 
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Mod 
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Mod 
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4539 

Mod 

0 

Pu-240 

94240 

4540 

Mod 

0 

* 
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GAMMA  PRODUCTION  DATA 
ALSO  AVAILABLE 

Yes 

No 

No 

No 

Yes 

Yes 

* 

Yes 

Yes 

No 

Yes 

Yes 

Yes 

No 

Yes 

Yes 

Yes 

Yes 

No 

Yes 

Yes 

No 

No 

No 

Yes 

Yes 

Yes 

Yes 

Yes 

Yes 

Yes 

No 

No 
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TABLE  II  - REGION  DESCRIPTION 
(region  numbers  are  circled  on  Figure  3.4c  for  clarity) 


* 


** 


REGION 

1 

2* 

3 

4** 

5 

6 

7 

8 
9 

10 


INPUT 

+1  -2 

+2  -3  -4  -5  -6  -7  -8  -9 
OR  +3  -5  -6  -10  OR  +4 
+5  -6 
+6 

+7  -9 
+8 

+9  -2  -7 
+7  +9 
+10 


Bodies  5 and  6 project  outside  body  3 in  the  X-Y  plane 
(see  Appendix  P) 

Bodies  5 and  6 buttress  each  other.  The  negation  of 
body  6 here  eliminates  possibility  of  machine  round 
off  error.  Negation  might  not  be  necessary  on  CDC 
machines. 
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One  possible  type  of  input  error  is  shown  in  figure  3.5b 
wherein  the  overlap  area  remains  undefined.  Another  type  of 
error  is  shown  in  figure  3.5c  wherein  the  overlap  area  has  been 
double-defined.  The  correctness  of  this  region  description  can 
be  tested  by  choosing  a single  point  in  the  region  and  testing 
the  region  description  at  this  point. 

The  checking  routines  provide  a general  method  of  obtaining 
a sufficient  number  of  check  points  within  a given  combinatorial 
geometry  description  and  the  checking  procedures  to  test  the 
single  definition  of  the  problem  space.  It  is  based  upon  the 
following  general  theorem  relating  the  minimum  number  (and 
location)  of  points  at  which  the  regions  are  singly  defined  to 
the  specification  of  the  entire  geometry  as  the  sun  of  singly 
defined  regions.  Consider  each  pair  or  triplet  of  intersecting 
bodies.  From  each  of  the  bodies  choose  a bounded  surface  such  that 
a corresponding  pair  or  triplet  of  intersecting  surfaces  occurs. 
Choose  single  points  infinitesimally  displaced  from  the  point  of 
intersection  within  each  of  the  volumes  enclosed  by  neighboring 
bounded  surfaces.  For  each  non-intersecting  body,  choose  two 
points,  one  immediately  interior,  one  immediately  exterior  to 
the  body.  If,  and  only  if,  all  these  points  lie  within  singly 
defined  regions  for  all  singlet,  pairs  and  triplets  of  bodies 
and  their  corresponding  bounded  surfaces,  the  problem  space  is 
uniquely  defined. 

In  the  SAM-CE  technique,  the  determination  of  intersection 
points  is  found  either  by  ray  tracing  (by  standard  methods) 
along  the  edges  formed  by  the  intersectin  of  two  body  planes* 

v 

*Ray  tracing  along  body  edges  is  also  utilized  by  the  code  to 
establish  whether  two  bodies  composed  entirely  of  planes 
buttress  or  truly  overlap  each  other. 


or  by  direct  solution  for  the  intersection  point  of  three  sur- 
faces. It  is  easily  established  that  for  two  body  intersections, 
ray  tracing  can  be  utilized  for  all  but  one  case  - a'  bounded 
quadratic  surface  of  one  body  intersecting  a bounded  surface  of 
the  second  body  in  a curve  that  does  not  abut  upon  an  edge  formed 
by  two  body  planes.  For  three  body  intersections  ray  tracing 
cannot  be  used  for  the  intersection  point  of  three  bounded  sur- 
faces at  least  two  of  which  are  quadratic. 

The  method  which  has  been  developed  and  coded  into  SAM-CE, 
solves  for  the  intersection  of  three  surfaces;  all  of  which  may 
be  quadratic.  The  procedure  reduces  the  problem  of  three  quad- 
ratic equations  into  one  which  invblves  finding  the  real  roots 
of  an  eight  degree  polynomial.  For  two  surface  intersections, 
a minimization  procedure  is  used  to  determine  a third  inter- 
secting surface  in  order  to  obtain  intersection  points.  Once 
the  intersection  points  are  found,  ray  tracing  methods  are  used 
to  determine  whether  they  actually  lie  upon  all  (two  or  three) 
bodies. 

In  the  interests  of  shortening  computer  running  time,  the  proced- 
ure utilizes  whenever  possible,  the  less  time  consuming  ray 
tracing  technique  as  opposed  to  the  surface  intersection  method. 

For  the  same  reason,  it  eliminates  from  consideration  all  body 
pairs  and  triplets  that  can  easily  be  shown  not  to  overlap. 

These  goals  are  obtained  by  the  following  sequence.  Place 
"Equivalent  RPPs"  (axes  parallel  to  the  coordinate  axes)  about 
each  body.  Test  these  for  two  body  overlap.  If  these  RPPs  do 
not  overlap  no  further  testing  is  necessary.  If  RPP  overlap 
does  occur  the  enclosed  bodies  themselves  may  overlap.  Place 


close  fitting  ARBs  or  BOXs  about  bodies  and  test  for  two  body 
overlap.  If  this  additional  check  still  indicates  possible 
overlap,  the  longer  procedure  is  required  and  the  code  tests 
for  true  twd  body  intersections  (by  ray  tracing  where  possible) . 
Test  triplet  intersections  for  only  those  bodies  which  intersect 
pairwise. 


Angular  Distribution 


Sources  may  be  either  isotropic  or  monodirectional  but  the 

same  angular  distribution  must  be  used  in  all  source  regions.*  (It 

should  be  noted,  however,  that  a source  may  be  generated  in  a 

finite  cone  by  specifying  an  isotropic  distribution  and  using 

exceedingly  high,  (e.g.,  1010)  angular  weights  to  kill  particles 

** 

which  are  generated  outside  the  desired  cone.)  At  the  present 
time,  the  source  must  be  isotropic  if  uncollided  flux-at-a-point 

estimates  (Section  3.2.8)  are  to  be  made. 

Energy  Distribution 

The  code  has  built  into  it  the  Cranberg  fission  neutron  spec- 


trum. If  this  option  is  selected,  no  energy  spectrum  input  is 
required.  If  an  arbitrary  spectrum  is  desired,  the  input  must 
contain  the  desired  energy  mesh  and  the  integrated  source  above 

each  energy  point  (i.e.,  a table  of  E vs  A maxS(E)dE  is  required). 
Monoonergetic  sources  may  also  be  specified. 

Time  Distribution (Time-dependent  problems  only) 

If  a time -dependent  problem  is  to  be  run,  the  user  must  supply 
a table  of  time  values  and  the  integrated  source  up  to  each  time 
(i.e.,  t vs  S(t)dt).  As  an  option,  all  radiation  may  be’  emitted 
’’t  a user-specified  single  value  of  time. 

3 . 2 . 4 . 2 External  Source 


The  user  may  supply  an  externally  generated  source  tape. 
Each  source  particle  must  be  described  by  14  parameters  in  the 
following  order: 


* Angular  biasing  of  an  isotropic  source  is  permitted,  however 


see  section  3.2.10  for  discussion  of  angular  weights 


3 Cartesian  Coordinates  (cm.)  (floating  point) 

3 direction  cosines  (floating  point) 

1 energy  (ev)  (floating  point) 

1 geometric  region  number  (integer) 

1 time  (sec.)  (floating  point) 

1 unused  parameter;  enter  a zero  here. 

1 weight  (floating  point) 

1 history  number  (integer) 

1 extra  carry  along  weight;  enter  a 1 here  (floating  point) 

1 particle  type  index;  enter  a 1 here  (integer) 

These  data  must  have  been  written  previously,  by  a separate  program, 
on  a file  designated  as  tape  15,  as  a series  of  binary  records  for  35 
particles  at  a time;  i.e.,  (14  parameters  per  particle)  x (35  parti- 
cles/record) =490  parameters  per  record. 


Given  the  array  A (14, 35)  the  proper  output  statement  is* 
WRITE (15)  ( (A (I , J) ,1-1,14) , J=1 , 35) . 

3. 2. 4. 3 Source  from  Previously  Generated  Interaction  Ta 


Using  a neutron  interaction  tape  (see  Section  3.2.1  ),  and 
gamma  ray  production  data  supplied  by  SAM-X,  SAM-F  can  generate, 
internally,  sources  of  secondary  gamma  radiation.  Note  that  the 
interaction  tape  was  generated  by  a previous  SAM-F  calculation  and 
has  a form  similar  to  that  of  the  external  source  tape  of  Section 
3. 2. 4. 2.  (See  Table  3.1,  on  page  207,  for  a complete  description.) 


* Experienced  programmers  prefer,  of  course,  WRITE (15)  A 
which  is  considerably  faster. 
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3.2.5  Output  Energy  Mesh 

During  tracking,  the  code  stores  fluxes  in  each  region  in  a 
set  of  energy  output  bins  specified  by  the  user.  The  number  (&100) 
and  width  of  these  bins  are  arbitrary.  The  bin  limits  must  be 
given  consecutively  in  the  input,  starting  with  the  highest  energy. 
The  upper  and  lower  bin  limits  for  all  supergroups  must  be  preceded 
by  minus  signs  as  explained  in  Section  3.2.3.  Care  should  be 
taken  to  insure  that  the  upper  energy  bin  limit  is  equal  to  or 
greater  than  the  highest  source  energy  to  be  generated  in  the  prob- 
lem. A cutoff  energy  is  also  specified,  which  instructs  the  code 
to  cease  tracking  any  particle  which  degrades  below  this  energy.* 

The  user  should  be  certain  that  the  lowest  energy  bin  limit  is 
lower  than  the  cutoff  energy.  In  essence,  there  must  be  a bin 
available  to  store  every  possible  energy  in  the  problem. 

(In  order  to  simplify  matters  for  the  user.  Appendix  E gives 
the  hierarchy  table  for  all  input  energy  limits  and  tables.) 

3.2.6  Time  Dependence 

SAM-F  enables  the  user  to  compute  particle  fluxes  as  a function 
of  time  as  well  as  energy  and  position.  The  u^er  selects  any  de- 
sired time  bin  structure  for  the  problem  and  enters  the  bin  limits 
in  consecutive  order  on  the  input  forms,  starting  with  the  highest 
bin.  Output  fluxes  will  be  given  in  this  bin  structure  in  the  edit. 

* However,  see  Section  3.2.11,  for  discussion  of  the  thermal 

neutron  option. 


147 


i 

t 

IV 

IT 

4 


•i 


The  longest  time,  used  in  specifying  the  time  bins,  will  be  used 
by  the  code  as  the  time-cutoff  parameter.  During  the  tracking 
process  the  code  computes  the  flight  time  of  a particle  between 
collision  points  from  its  velocity  (or  energy) . All  nuclear  inter- 
actions are  assumed  to  occur  instantaneously.  By  accumulating  the 
flight  times  for  each  particle,  the  code  is  capable  of  storing 
particle  fluxes  in  the  proper  output  time  bins. 

3.2.7  Scoring  Regions 

A scoring  region  is  one  in  which  a flux  contribution  is  com- 
puted for  each  particle  which  passes  through  it.  In  a nonscoring 
region  no  such  computation  is  made,  so  that  the  output  edit  provides 
fluxes  only  in  those  regions  designated  in  the  input  as  scoring 
regions . 

In  some  problems  it  is  desired  to  know  the  flux  in  every  region 
separately,  in  which  case  each  region  in  the  problem  would  be  de- 
fined as  a scoring  region  with  a different  number.  In  some  prob- 
lems, however,  two  or  more  regions  may  be  completely  symmetric 
with  respect  to  the  source,  in  which  case  the  fluxes  in  these  sym- 
retric  regions  could  be  combined  without  any  loss  of  information, 
and  in  fact,  an  improvement  in  the  accuracy  will  be  obtained.  Each 
of  these  regions  then  would  be  designated  by  the  same  scoring 
region  number.  In  still  other  problems  it  may  be  unnecessary  to 
know  the  fluxes  in  certain  regions.  These  should  then  be  given 
scoring  region  number  zero,  which  tabs  them  as  nonscoring. 
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FLOW  CHART  (Simplified) 


Overlay  AGEOM 


option 


GCHRPP 

(checks  for  un- 
defined and  multiply 
defined  regions) 


FLOW  CHART  (Simplified) 
Overlay  ADATA 


3.3.1  SAM-F  Routines  (Alphabetical)  - General  Description 
This  section  gives  a brief  description  of  every  routine  in 
the  program,  given  in  alphabetical  order.*  These  brief  descrip- 
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tions  are  followed,  in  Section  3.3.2,  by  detailed  descriptions 
of  those  important  routines  which  are  flagged  by  an  asterisk  (*) 


superscript. 

Routine 

ADATA 


AGEOM 


AMONTE* 


AMOUT 


AN I SOT 


ARG* 
ARP REP 


ASORTT* 


Description 

The  controlling  program  of  the  second 
primary  overlay.  It  oversees  three 
secondary  overlays  which  read  in  the 
input  data. 

The  controlling  program  of  the  firBt 
primary  overlay.  In  this  overlay, 
geometry  data  is  read  in,  organized, 
and,  as  an  option,  checked  for  logical 
errors. 

The  controlling  program  of  the  fourth 
primary  overlay.  It  oversees  the  routines 
which  generate  source  particles,  do  the 
tracking  and  scoring  and  perform  the 
tallying  of  results.  AMONTE  also  controls 
the  shuttling  in  and  out  of  banded  cross 
section  data  and  the  organization  of  the 
particle  histories  into  aggregates. 

The  controlling  program  of  the  fifth 
primary  overlay.  It  calls  SUBED  to 
perform  the  final  edit. 

Called  by  DR3  to  determine  the  center  of 
mass  angle  of  scattering  for  anisotropic 
events . 

The  general  angle  reseleOtion  routine  of 
the  bounded  f lux-at-a-point  procedure. 

Called  by  ARG,  as  an  aid  in  angle  reselect- 
"ion  for  two  special  cases:  (1)  neutron 
energy  fust  above  the  inelastic  threshold 
or  (2)  for  neutron  scattering  by  hydrogen. 

The  third  primary  overlay.  For  secondary 
gamma  ray  problems,  this  routine  presorts 
the  precursor  interaction  events. 


AUTOBA  Calculates  and  supplies  to  BAND  the  energy 

limits  for  cross  section  bands. 

*A11  routines  involved  in  the  geometry  package  begin  with  the  letter 
"G".  All  other  routines  begin  with  letters  other  than  "G". 
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FILL 

FLUP* 

FLUPV* 

GARB 

GARBER 

GARBIH 

GARBOX 


The  edit  routine.  It  prints  out  flux 
and  flux-dependent  (e-g.,  dose)  answers 
as  functions  of  energy,  time  and  scoring 
regions.  The  routine  also  prints  out 
"point"  and  "small  volume"  detector  results. 

An  auxiliary  routine,  called  by  BAND,  used 
to  fill  in  the  cross  section  arrays  and  the 
scattering  data  tables  in  the  banded  arrays. 

The  flux-at-a-point  estimation  routine. 

The  routine  calculates  scores  as  functions 
of  energy,  time  and  detector  point.  The 
scores  are  stored  in  the  MASTER  array  and 
are  printed  out  by  EDIT. 

Similar  to  FLUP.  Used  for  "small  volume" 
detectors. 

Special  geometry  routine  to  process  the 
input  for  the  ARB  body  of  the  Combinatorial 
Geometry  package. 

Geometry  checker  routine  called  by  GARBOX 
to  store  planar  surface  parameters  of  RPP, 
BOX,  RAW,  in  the  MASTER  array  for  later  use. 

Geometry  checker  routine.  Called  by  GCHECK 
to  convert  ARB  plane  parameters  into  ARB 
point  parameters. 

Geometry  checker  routine  called  by  GCHRPP 
to  place  closest  fitting  BOX  or  ARB  about 
each  body  containing  a quadratic  surface. 

It  stores  these  data,  the  normal  vector, 
and  distance  from  origin  to  the  planes  of 
a body  containing  a quadratic  surface  for 
use  in  region  checking.  GARBER  is  called 
to  store  the  plane  parameters  for  planar 
bodies . 


GATHER 


GCHECK 


GCHKBD 


Geometry  checker  routine  called  by  GTSURF 
and  GTQDPT  to  gather  terms  of  a rotated 
and  translated  "principle  axis"  equation 
to  form  coefficients  of  the  general  trans- 
formed quadratic  surfaces. 

Main  geometry  checker  routine  called  by 
AGEOM  for  routines  (GCHKBD,  GCHKRG,  GCHRPP 
and  GEXTRA  - see  below)  to  perform 
geometry  checking  of  various  tynes.  It 
also  calls  GARBIN  (see  above)  to  store 
point  data  for  ARBs  needed  for  checking. 

Geometry  checker  routine.  Performs  various 
body  data  checks.  Calls  GEQUIV  (see  below) . 


GCHKRG 


GCHRPP 


Geometry  checker  routine  called  by  GCHECK 
to  check  for  simple  errors  in  region 
descriptions  and  to  provide  a table  of 
region  references  for  each  bodv. 

Geometry  checker  routine.  Once  the  body 
and  simple  region  data  have  been  checked 
by  GCHKBD  and  GCHKRG,  GCHRPP  examines  the 
RPPs  placed  about  each  body  to  establish 
possible  two  body  overlaps.  If  RPP  over- 
lapping occurs,  GARBOX  is  called  to  place 
an  equivalent  BOX  or  ARB  about  each  non- 
planar  body  of  the  pair  and  to  store  these 
data  and  body  plane  parameters  in  the 
PIASTER  array.  Overlapping  of  the  equivalent 
body  pairs  and  the  two  body  intersections 
(involving  a planar  body)  are  then 
checked  (via  ray  tracing  along  body  edges) 
by  a call  to  GIRTWB.  For  overlapping  body 
pairs,  a call  to  GIRTWC  then  checks  the 
remaining  two  body  intersection  points  for 
no  or  multiple  region  specification.  For 
intersecting  pairs,  GARBOX  is  called  for 
third  bodies  which  intersect  each  of  the 
bodies  of  the  pair.  GIRTWC  then  checks 
three  body  intersections  for  improper 
region  specification.  Printout  of  the  two 
and  three  body  overlaps  (TBO) , buttressing 
surfaces  for  two  planar  bodies  (BSTB) , 
equivalent- BOX  or  ARB  overlap  (EBOXARBO) , 
equivalent  RPP  overlap  (ERRPO)  occurs  for 
values  of  IORPP=0...3  respectively.  Higher 
values  of  IORPP  lead  to  printout  of  the  new 
case  in  addition  to  those  for  lower  numbers. 


GCROSS 


GDETER 


Used  by  checker  routines  to  compute  cross 
products . 

Geometry  checker  routine  called  by  GTPLPT 
to  calculate  the  determinant  of  a 3x3 
matrix. 


GDOT 


GENI 


Used  by  checker  routines  to  compute  dot 
products . 

The  major  geometry  input  processing  routine 
called  by  AGEOM.  The  routine  first  reads 
body  geometry  data  then  calls  GFREE  to  read 
the  region  description  (in  a free  field 
format)  and  puts  the  data  into  the  FPD  and 
MA  arrays  in  the  form  required  by  the 
tracking  routines. 
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GEQUIV 

GETIR 

GETXB 


GEXTRA 


GFREE 


GG 


GGTEC 


GIROB 


CIRTWB 


Used  by  checker  routines  to  place  an  RPP 
about  each  body  in  the  geometry. 

Used  to  locate  the  region  number  of  a 
given  spatial  location.  During  geometry 
checking  it  identifies  points  in  undefined 
and  multiply-defined  regions  (see  Fig.  3.5). 

Geometry  checker  routine  called  by  GTQDPr 
to  test  if  a point  lies  on  each  of  two  or 
three  bodies.  It  first  calls  GTBDPT  to 
obtain  an  interior  point  for  each  of  the 
bodies  and  then  tests  to  see  if  the  distance 
from  this  point  to  the  body  surface  is  the 
same  as  distance  from  this  point  to  the 
point  of  interest.  If  so,  GTIRNT  is  called 
to  check  point  for  region  correctness. 

Geometry  checking  routine.  Calls  GPOINT 
(see  below)  and  also  sets  up  the  two- 
dimensional  point  grid  for  picture  printouts. 

This  routine  is  called  by  GENI  to  read 
region  data  in  a free  field  format. 

The  distance  calculating  routine.  Given  a 
position  X and  direction  W and  a body  number, 
this  routine  computes  the  two  distances  RIN, 
ROUT  measured  from  X to  the  body. 

Special  distance  calculating  routine  for  the 
truncated  elliptical  cone  (TEC) . 

Geometry  checker  routine  called  by  GCHRPP. 
Used  to  locate  points  to  check  for  possible 
undefined  and  multiply-defined  regions  near 
the  surface  of  individual  bodies. 

Geometry  checker  routine.  Called  by  GCHRPP 
to  find  intersections  of  two  bodies  by 
tracing  along  or  parallel  to  the  edge 
formed'  by  two  planes  of  a planar  body.  It 
operates  in  two  modes.  In  the  first,  ray 
tracing  along  the  edges  of  an  equivalent 
ARB  or  BOX  of  a body  with  a quadratic  surface 
determines  any  overlap  with  a second  equival- 
ent ARB  or  BOX.  In  the  second  mode,  the 
intersection  points  of  a planar  body  with 
another  body  (not  necessarily  planar)  are 
obtained  by  ray  tracing  in  the  four  quadrants 
formed  by  two  intersecting  planes.  By  this 
means  intersection  points  and  buttress  points 
(buttress  points  are  obtained  for  two  planar 
bodies)  are  obtained "for  the  two  bodies. 
Identification  of  undefined  and  multiply- 
defined  region  specification  at  these  points 
is  obtained  by  calls  to  GETIR. 
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GIRTWC 


Geometry  checker  routine.  Called  by  GCHR^P 
to  check  two  and  three  body  intersection 
points  which  are  found  by  rav  tracing  (by  a 
call  to  GG)  along  all  edges  formed  by  all 
possible  plane  surfaces  taken  one  from  each 
of  two  of  the  bodies.  It  uses  (at  a starting 
point  outside  the  bodies)  four  rays  parallel 
to  the  edge,  one  in  each  of  the  four  Quad- 
rants formed  by  the  two  planes,  in  order  not 
to  miss  buttressing  surfaces.  If  a point 
along  a ray  falls  on  all  the  bodies,  a call 
to  GTIRNT  checks  the  point  for  correct 
region  specification.  For  two  body  checking 
the  check  point  chosen  along  the  ray  is  at 
the  entry  of  the  ray  into  the  second  body 
after  passage  into  the  first  body  surface. 

For  three  bodies,  ray  tracing  is  done  along 
edges  formed  by  planes  of  two  of  the  bodies. 
Therefore,  both  the  last  entry  and  first 
exit  points  along  the  ray  are  region  checked. 
By  this  means,  points  along  all  possible 
two  plane  surface  intersections  are  obtained. 
Further  checking  occurs  via  a call  to  GTQDPT. 


An  auxiliary  routine,  called  as  an  option 
by  Gl,  for  debug  printout. 

The  main  geometry  tracking_routine.  Given 
a position  X,  a direction  W,  and  a region 
IR,  the  routine  will,  calculate  the  distance 
"S"  from  the  point  X to  the  next  region  in 
the  direction  V7.  The  routine  also  determines 
IR',  the  next  region  to  be  encountered. 


GPACK 


GPOINT 


GRAFIC 


GRID 


Geometry  checker  routine.  Called  by  GCHRPP 
to  pack  overlap  and  intersection  information 
of  body  pairs  in  the  MASTKR  array. 

Geometry  checker  routine.  Locates  and  prints 
out  the  region  of  all  user-supplied  points. 

Prints  out  the  two-dirtensional  pictures. 

Finds  the  proper  region  number  for  all 
points  in  the  two-dimensional  grid  for 
picture  display. 


GTARBX 


Geometry  checker  routine.  Called  by 
GIRTV7B  to  store  equivalent  ARB  or  BOX  for 
bodies  containing  a quadratic  surface  in 
unused  locations  of  the  FPD  and  MA  arrays. 


GTBDPT 


Geometry  checker  routine.  Called  by  GETXB 
to  obtain  interior  point  of  a body. 


174 


I 


I 

GTIRNT 

► '■ 

5 A 

V J 

-1 

| . GTPLPT 

[•  : 

[ j 

GTQDPT 

- 

i 

1 

‘I 


GTRQA5 

GTRQB3 

GTRQFN 

GTRQLN 

GTRQMX 

GTRQPL 

GTRQR3 


GTRQSF 

I 


Geometry  checker  routine.  Called  by  Girtwc 
and  GTQDPT  to  set  RIN  and  ROUT  parameters 
stored  for  each  of  the  two  or  three  inter- 
secting bodies  so  that  a subsequent  call 
to  GETIR  effectively  involves  points  in 
each  of  the  quadrants  or  octants  formed 
by  two  or  three  intersecting  surfaces. 

Geometry  checker  routine.  Called  by  GIRTWC 
to  obtain  a space  point  on  the  ed$e  formed 
by  two  planes. 

Geometry  checker  routine.  Called  by  GIRTWC 
to  obtain  intersection  points  of  all  pairs 
or  triplets  of  surfaces  from  each  of  two  or 
three  bodies,  no  more  than  one  of  which  is 
a plane,  via  calls  to  GTRQU . All  possible 
intersection  points  are  found  for  three  sur- 
faces and  a minimization  procedure  called 
by  GRTQU  is  used  to  obtain  a sufficient  number 
of  check  points  prior  to  calling  GTRQU.  GTSURF 
and/or  GATHER  are  called  to  obtain  surface 
parameters  from  the  stored  body  parameters  for 
two  surfaces.  A call  to  GETXB  determines 
whether  these  points  lie  on  all  the  bodies. 

If  so,  GETXB  calls  GTIRNT  to  check  the 
correctness  of  the  region  descriptions. 

Called  by  subroutine  GRTQU  to  perform  a 
partial  reduction  of  three  quadratic 
equations  of  three  unknowns. 

CAlled  by  subroutine  GTRQU  to  perform  a 
partial  reduction  of  three  quadratic 
equations  of  three  unknowns. 

Called  by  subroutine  GTRQPL  to  evaluate  a 
polynomial. 

Called  by  subroutine  GTRQU  to  determine 
whether  the  equations  being  solved  are 
linearly  independent. 

Called  by  subroutine  GTRQU  to  perform  a 
random  four-dimensional  rotation. 

Called  by  subroutine  GTRQU  to  obtain  the 
roots  of  a polynomial. 

Called  by  subroutine  GTRQU  to  perform  the 
final  reduction  of  three  quadratic  equations 
of  three  unknowns  to  a single  eight  degree 
equation  of  a single  unknown. 

Called  by  subroutine  GTRQU  to  set  up  a 
random  four-dimensional  rotation  matrix. 
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GTRQU 


GTSURF 


Called  by  GTQDPT  to  determine  the  coordinates 
of  the  intersection  of  three  quadratic  sur- 
faces, or  of  selected  points  on  the  inter- 
section of  two  quadratic  surfaces. 

Geometry  checker  routine.  Called  by  GTQDPT 
to  obtain  surface  parameters  for  bodies 
having  quadratic  surfaces.  The  method  used 
is  to  consider  all  terms  of  a quadratic 
surface  with  a specific  position  and  orien- 
tation. The  contribution  of  these  terms  to 
the  quadratic  of  interest  is  obtained 
utilizing  rotation  and  translation  of  coordin- 
ates via  a call  to  GATHER. 


GUN  IT 


Used  by  checker  routines  to  compute  unit 
vector. 


GVECTR 


INPUTD 


IMPUTE 


LENTHS 


MAIN* (SAMF) 


NDQSET 


PICK* 


RANDV7 


Geometry  checker  routine.  Checks  angles 
made  by  various  vectors  in  the  description 
of  the  geometric  bodies. 

Reads  in  most  of  the  Monte  Carlo  input  data 
and  stores  it  in  the  MASTER  array  for  use 
by  the  calculation  routines.  Region  speci- 
fication, importance  sampling  data,  output 
energy  meshes,  and  detector  data  are 
handled  by  this  routine. 

An  auxiliary  routine  used  to  flag  very 
distinctly  all  errors  in  the  Monte  Carlo 
input  data. 

An  auxiliary  routine  called  by  BAND  and 
FILL.  Used  to  determine  lengths  of  energy 
tables  in  a given  energy  hand. 

The  overall  control  routine  for  SAM-F. 

The  size  of  the  MASTER  array  is  set  by  SAMF. 

An  auxiliary  function  to  find  the  length  of 
specified  arrays. 

The  control  routine  for  superbin  la tents. 
This  routine  stores  particles  whose  energies 
are  not  included  in  the  superbin  currently 
being  processed.  These  particles  are  stored 
by  PICK  in  either  the  central  memory  or  on 
tapes  depending  on  the  amount  of  core  avail- 
able for  the  particular  problem. 

An  auxiliary  routine  called  by  FILL  to  read 
binary  blocks  of  data  from  the  element  data 
tape  (EDT)  and  to  write  binary  blocks  of 
reduced  data  on  the  banded  organized 'data 
tape  (ODT) . 
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RANF 


A built-in  random  number  generator. 


REED 

SEEK* 

SOUCAL* 

SOUGAM* 

SOUGEH* 

SOUPIC* 

SOUSEC 

SUBED 

TALLY* 

TDIFIN 


An  auxiliary  routine  to  road  in  binary 
arrays,  or  partial  arrays,  very  rapidly. 

Given  an  array,  A,  v/ith  elements  monotoni- 
cally  decreasing  and  a variable,  X,  this 
routine  will  search  through  the  A array 
and  determine  the  bin  containing  X. 

The  input  processor  for  the  source  infor- 
mation required  by  the  Monte  Carlo  routines 
for  primary  transport  problems.  The  routine 
processes  energy  and  time  spectra  data,  and 
source  region  data  which  are  stored  in  the 
MASTER  array. 

The  routine  which  generates  the  source 
particles  for  secondary  transport  problems. 

The  routine  which  generates  the  source 
particles  for  primary  transport  problems  - 
except  for  the  external  source  option. 

The  routine  which  calls  either  SOUGAM  or 
SOUGEN  to  generate  secondary  or  primary 
source  particles,  respectively.  It  will 
also  read  in  source  data  from  an  external 
source  tape.  In  addition,  SOUPIC  controls 
the  angle  of  emission  reselection  procedure 
for  problems  involving  one  or  more  point 
detectors. 

The  input  processor  for  the  gamma  product- 
ion data  for  secondary  transport  problems. 
Creates  file  13,  the  organized  gamma 
production  data  tape. 

The  main  edit  routine.  The  routine  reads 
the  answer  arrays  from  the  statistical 
aggregate  tape  and  prepares  it  for  editing 
by  the  EDIT  routine.  Flux-dependent 
responses  are  also  calculated  in  this 
program. 


> 


A summary  routine.  The  routine  prints  a 
one  line  summary  of  results  for  each 
statistical  aggregate.  Quantities  such 
as  number  of  collisions,  absorptions,  and 
energy  deposition  are  printed  for  each 
aggregate . 

Celled  by  INPUT!)  if  thermal  neutron  diffu- 
sion option  has  been  involved.  TDTFI**  reads  ) 

and  preprocesses  the  diffusion  option  input. 


TERP 

THERM 

( 

V„L  ^ 

TIMEX 

TRALA* 

TRALAV* 

TROPIC 

VCALC 

WRIT 

WRT14 


A linear  interpolation  routine. 

A subroutine,  called  bt  AUTOBA,  which  pre- 
calculates cross  sections  avcraqed  over 
the  Maxwellian  spectrum  for  the  thermal 
neutron  option. 

Used  to  keep  track  of  computer  running  time 
during  Monte  Carlo  execution. 

A tracking  routine  usdd  by  FLUP,  the  main 
f lux-at-a-point  routine.  It  tracks  from  an 
initial  collision  (or  source)  point  to  a 
detector  point. 

Similar  to  TRALA.  Used  by  FLUPV,  the  flux- 
in-a-small-volune  routine. 

A routine  to  generate  a vector  of  direction 
cosines  from  an  isotropic  distribution. 

The  volume  computation  routine.  Region 
volumes  are  computed  by  numerical 
integration. 

An  auxiliary  routine  to  write  out  binary 
arrays,  or  partial  arrays,  very  rapidly. 

A routine  to  write  14-word  records  onto 
tape.  The  routine  is  called  whenever  a 
transmission  or  interaction  is  to  he  put 
on  tape. 
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Subsequent  to  the  direction  checking  stage,  at  most  one  live 
detector  sphere  is  intersected  by  the  originally  sclectod  rav.  if 
such  an  intersection  has  survived,  a new  direction  is  selected  in 
the  reselection  staae. 

The  final  stage  in  the  reselection  procedure  involves  the  cal- 
culation of  a weight  adjustment  factor  which  compensates  for  the 
biasing  introduced  by  reselection. 

Program  ASORTT 

ASORTT  reorganizes  the  neutron  interaction  tape  so  as  to  pro- 
vide (for  secondary  problems)  a memory  saving  scheme  analogous  to 
the  BAND  feature  for  transport  cross  sections.  Whereas  BAND  auto- 
matically splits  cross  section  data  into  two  or  more  energy  ranges, 
to  be  treated  one  at  a time,  ASORTT  makes  provisions  so  that  during 
the  generation  of  secondary  gamma  rays,  gamma  production  data  for 
one  element  only  resides  in  the  MASTER  array  at  any  given  time. 

The  basis  for  this  memory-saving  process  lies  in  the  prior 
creation  of  an  organized  gamma  production 'data  tape  (OGPDT)  by 
SOUSEC  and  a sorted  interaction  tape,  (SIT) , by  ASORTT. 

The  format  of  the  OGPDT  is  similary  to  the  original  gamma 
production  data  tape  described  in  Appendix  D,  except  that: 

a)  the  first  two  words  for  each  element  are  dropped, 

b)  the  pointers  are  made  relative  to  the  starting  address 
in  the  PIASTER  array, 

c)  the  data  for  each  element  are  written  as  one  logical 
record  on  the  OGPDT. 
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AMONTE,  the  controlling  Monte  Carlo  overlay  program  first 
considers  a potential  new  ray.  This  can  be  either  a source  parti- 
cle (read  from  source  tape  or  generated  internally  by  subroutine 
SOUPIC)  or  a picked  up  latent  from  a previous  interaction.  AMONTE 
next  tests  the  energy  of  the  considered  particle.  If  it  is  below 
the  lower  bound  of  the  superbin  currently  treated,  it  is  stored 
as  a latent.  If  the  energy  of  the  particle  is  within  the  current 
superbin,  it  is  transmitted  to  CARLO. 

CARLO  proceeds  to  track  the  particle  and  to  score  answers. 
Collisions  made  by  the  particle  are  stored  by  CARSCA.  The  track- 
ing procedure  in  CARLO  is  as  follows. 

Given  the  region  number,  IR,  the  energy,  E,  and  the  direction 
of  flight,  WB,  a sampling  weight* 


w = WlR  x WE  x Wfi 

is  calculated.  A test  is  made  to  determine  whether  or  not  the 
thermal  diffusion  option  is  being  invoked. 

If  the  thermal  diffusion  option  is  not  being  invoked,  the 
"ordinary"  tracking  procedure  is  used.  First,  subroutine  DR1  is 
called,  which  provides  the  total  macroscopic  cross  section  in  the 
region  IR. 

When  point  detectors  are  being  used,  (and  if  the  extended 
particle  path  passes  sufficiently  close  to  a specified  detector) 
preliminary  geometrical  calculations,  necessary  to  bias  the 
collision  positions,  are  carried  out.  This  biasing  is  that 
required  to  make  flux  estimates  for  point  detectors  bounded. 


* See  Appendix  F for  a detailed  discussion  of 

the  use  of  importance  sampling  in  Monte  Carlo  calculations. 


182 


The  subroutine  is  called  when  a particle  enters  collision  at 
some  point  within  a region  of  given  composition.  Its  first 
function  is  to  select  the  nuclide  in  the  composition  with  which 
tne  particle  collides.  The  probability  that  a particle  with  energy 
E will  intereact  with  nuclide  k of  the  composition  R is  given  by 

p(k,R)  = C(kfR)oT(k,E)/yT>R(E)  f p(K#R)  « ? p(n,R) 

n-1 

where  y_  -(E)  is  the  total  macroscopic  cross  section  of  the  compo - 

1 # K 

sition,  and  C(k,R)  is  the  concentration  of  nuclide  k in  region  R. 

Pk  -24 

C(k,R)  = N x 10 

"k  ° 

3 

where  pk  = density  of  the  nuclide,  k,  in  the  composition  (gm/cm  ) 

^ = its  atomic  weight 

24 

Nq  = 0.6023  x 10  ; Avogadro's  number 


By  normal  Monte  Carlo  procedure,  DR3  generates  a random 
number,  £,  and  successively  compares  £ with  P(k,R),  for  all 
values  of  k,  until  the  condition 
£<P (k,  R) 


is  satisfied.  Then  the  interaction  is  with  that  k-th  nuclide.  | 

Next,  DR3  determines  the  type  of  interaction  suffered  by  | 

the  particle.  The  probability  that  a particle  with  energy  E ! 

will  suffer  the  "i-th"  type  interaction  with  the  already  selected 
"k-th"  nuclide,  is  given  by: 

P(i,h)  » OifEj/o^E) 
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J 


path  from  S to, I.  If  the  ray  misses 

the  detector  completely  no  score  is  made  and  control 

is  returned  to  the  calling  program. 

If  interception  is  made  TRALAV  next  calculates  IJ 
the  flight  distance  (cm)  through  the  detector.  Control 
is  returned  to  FLUPV  where  a track  length  score  of 


WT  x g x e"A  x ij  x 


[41 


for  >d>0  or 


WT  x gg  x e x ij 


for  Ad=0  (vacuum) 


is  made,  where  *d  is  the  mean  free  path  flight  path 
through  the  detector  and  WT  is  the  weight  of  the 
particle  at  point  S,  prior  to  collision. 

Note  that  if  the  ray  starts  inside  the  sphere  it  is 
allowed  to  scatter  in  an  unbiased  direction  (non- 
source) or  to  proceed  in  its  original  direction 
(source)  and  the  same  estimate  as  above  is  made 
with  g =1,  and,  if  the  particle  is  inside  the 
detector  as  well  as  inside  the  sphere,  e”Vl.  if 
the  particle  is  inside  the  sphere  but  outside  the 
detector,  and  misses  the  detector,  no  score  is  made. 


The  mathematical  technique  to  select  the  scattering  angle  at 
S is  of  some  interest  and  is  now  presented. 
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where  3 is  a random  number  uniformly  distributed  between  0.  and  1. 

Once  cosa  has  been  selected  it  is  simple  by  usual  Monte  Carlo 
procedures  to  select  a random  azimuthal  angle  and  to  calculate  the 
direction  cosines  of  SI*  which  are  designated  Wx2,  Wy2  an<*  WZ2* 

Then,  if  0 is  the  selected  angle  of  scattering: 

COS0  = Wx1*Wx2+Wy1*Wy2+Wz1*Wz2 

and 

gg  = g(cose)  x 2v  (1-cosi|i) 
where  g(cos0)=  for  isotropic  scattering. 

S inside  sphere,  but  outside  detector 

Select  WX2,  Wy2 , Wz2  from  an  isotropic  distribution.  Then 
gg=l  and  score  only  if  path  intersects  detector. 

S inside  detector 

Again,  select  WX2,  Wy2  and  WZ2  from  an  isotropic  distribution 
and  then  gs=l. 

Note:  In  order  to  increase  calculational  efficiency  (by  decreasing 

the  number  of  rays  which  intercept  the  sphere  but  miss  the  detector) 
MAGI  investigated  the  effects  of  splitting  the  cone  into  an  inner 
and  outer  cone.  The  inner  cone  goes  from  angle  0 to  <//2  and  the 
outer  cone  from  \\i/2  to  <ii , (see  above  diagram)  . The  code  fo.'ces  3 
particles  into  the  inner  cone  for  each  particle  that  goes  into  the 
outer  cone  - with  corresponding  weight  adjustments.  This  modifica- 
tion does  seem  to  increase  the  efficiency  of  the  technique  for 
cubic  and  cylindrical  detectors  (with  the  diameter i about  equal  to 
the  altitude) , and  is  currently  coded  in  SAM-F.  Investigations 
along  these  lines  should  continue. 
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Description  of  14-Wbrd  Particle  Descriptive  Array 


WORD(s)  NAME 


Stored  in  Latent 
Table  (Supergroup 
Option  Only) 


Stored  on  External 
Source  Tape 


Interactions  Stored 
on  Interaction/ 
Transmission  Tape 


Transnissi 
Stored  on 
Interaction/ 
Transmission 


Cartesian 

Coordinates 


Cartesian 

Coordinates 


Cartesian 

Coordinates 


Direction  of 
Particle  fran  the 
Source  or  Caning 
out  of  Collision 


Direction  of 
Particle  frcm 
the  source 


Energy  of  Particle  Energy  of  Particle 


Direction  of 
Precursor  Primary 
Particle 


Region  Number 


Region  Number 


Region  Number 
(Must  Match 
Precursor  Run) 


Region  Number 
(Must  Match 
Precursor  Run) 


=100x  KDLIV+IDLIV 
if  J12345=l 


=200000xKDLIV+SIGN 
(KDLIV)xIATWr 
Note:  = IAXWT 
for  no  detectors 


Statistical  Weight  Statistical  Weight  Statistical  Weight  Statistical  Wei' 


History  Nutber  History  Number 


Precursor  Primary 
Particle  History 
NUnber 


Carry-Along 

Weight 


Carry-Along 

Weight 


distance  to 
detector 
if  J12345=5 


is  ignored 
if  J12345=4 


1 for  source 
_ garticle_ 


2 for  particle 
caning  out  of 
elastic 
scattering 


3 for  particle 
coming  out  of 
inelastic 
scattering 


? for  small-uol 
latent 


5 for 


nt  3^ 
latent 


5 for  point  de- 
tector  estimal 
transmission 


3.3.3  Overlay  Structure  of  SAM-F 


The  following  is  the  overlay  structure  of  SAM-F. 


SAMF 

GETIR 

GG 

GGTEC 

GP 

G1 

NDQSET 

RANF 

REED 

SEEK 

TIMEX 

TROPIC 

WRIT 


AGEOM 


GCHECK 

GCHKBD 

GCHKRG 

GCHRPP 


GEXTRA 

GFREE 

GPOINT 

GRAFIC 


(All  otht 
"G"  rou- 
tines art 
controllt  d 
by  GCHRPI 
for  re- 
gion 

checking 


ADATA 

DSPLAY 

THERM 

VCALC 


DRGE 

DR1 

DR3 

FLOP 

FLUPV 

PICK 

SOUGAM 

SOTGEN 

SOUP I C 

TALLY 

TERP 

TRALA 

TRALAV 

WRT14 
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3.3.4  Glossary  of  Important  Common  Arrays 


The  following  is  a description  of  each  variable  in  blank 
COMMON  and  in  some  of  the  more  important  labeled  common  arrays. 


EOUT(IOO) 


EWTAB ( 50 ) 
ANGLE (50) 


"Blank"  Common 

An  array  containing  the  output  energy  bins  for 
flux  results.  The  array  contains  all  bins  for 
all  output  supergroups. 

The  energy  mesh  for  energy  importance  sampling. 
The  cosine  mesh  for  angular  importance  sampling. 


TTAB (50) 
MASTER (7, 000) 

KDUM 


The  time  mesh  for  iifie  dependent  problems. 

The  master  Storage  array  containing  input  and 
flux  data.  A complete  description  appears  in 
Appendix  H.  The  array  is  equivalenced  with 
array  ASTER. 

Dummy  variable  needed  by  function  MDQSFT  to 
change  the  memory  allocation  for  the  MASTER 
array . 


ISC 

NREG 

IRW 

IEW 

IAM 

IANG 


COMMON  REGPAR 

The  scoring  region  number  for  region  IR. 

The  composition  number  for  region  IR. 

The  location,  in  MASTER,  of  the  energy  weighj 
for  region  IR. 


tOrfT  in  MASTER,  of  the  energy  weight 
table  for  region  IR. 

The  location,  in  MASTER,  of  the  aiming  angle 
for  region  IR. 

The  location,  in  MASTER,  of  the  angular  weight 
table  for  region  IR. 
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2 . 4 . 1 . 1 Geometry 

r Input 

Item  1 General  Information  (Format  515, E15. 5, 10A4 

IGOPT  = 0,  Suppress  both  geometry  checking  and  picture 
taking  options 


= 1,  Check  for  body  and  region  description  errors 

= 2,  Make  picture (s)  of  plane  slice (s)  through  the 
geometry  and/or  check  user- specified  points 

=3,  Do  both  1 and  2 

* 

I SEC  Escape  Region.  (See  also  footnote  2.) 

* 

IORPP  = 0,  During  geometry  checking,  suppress  all  printout 
except  for  region  error  messages. 

=1,  During  geometry  checking,  provide  buttressing 

surface  (BSTB)  and  two  and  three  body  intersection 
(TBI)  information,  in  addition  to  printout  of 
IORPP=0 . 

= 2,  During  geometry  checking,  provide  equivalent  ARB 
or  BOX  overlap  (EBOXARBO)  information,  in  addi- 
tion to  printout  of  IORPP*l 

= 3,  During  geometry  checking,  provide  equivalent  RPP 
overlap  (ERPPO)  information L in__addit4on -to — 

printout'  of  IORPP=2 

IPRINT  = 0,  Print  out  body  and  region  data  which  follow 

= 1,  Print  out  body  and  region  data,  the  internal 
arrays  in  which  they  are  stored,  as  well  as 
additional  information  regarding  ARB  bodies. 

■ 2,  Suppress  all  geometry  printout 


These  items  need  be  entered  only  if  IGOPT-1  or  3.  Otherwise,  they 
may  be  left  blank. 

**The  escape  region  is  a special  region  which  generally  encloses  all 
other  regions.  In  ray-tracing  part  of  the  code,  rays  entering  the 
escape  region  are  terminated. 
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SCALE 


= 0,  no  printout 
= 1,  main  debug  printout  option 

- 2,  same  as  IDBG=1  plus  much  extra  geometry 

tracking  printout. 

- 3,  same  as  IDBG»1  plus  gamma  ray  production  data 

printout. 

■ 4,  same  as  IDBG»2  plus  gamma  ray  production  data 
printout. 

Scale  Factor.  Multiply  all  Combinatorial  Geometry 
dimensions  by  this  factor.  (Default«1.0) 

Title,  40  arbitrary  Hollerith  characters. 
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Item  3 Region  Cards  (Format  2X, A3, free  field  for  rest  of  card) 

Each  region  must  be  numbered  and  described  by  a logical  combin- 
ation of  the  bodies  which  make  up  that  region.  Use  as  many  cards 
as  necessary  to  describe  each  region  and  begin  each  region  on  a now 

Input 

Blank 

Arbitrary  non-blank  Hollerith  Data 
Free  field  input  describing  the  logical  combination 
of  bodies  making  up  the  region.  All  bodies  specified 
are  usually  preceded  by  either  the  "+"  or 
operators.  If  no  operator  is  given,  the  code  will 
assume  a "+".  (For  IBM  [029]  punched  cards  the  "+" 
must  not  be  explicitly  given.)  Insert  OP.  operators 

a 

if  needed*,  before  the  "+",  or  blank  operators*. 
Blanks  may  be  freely  left  anywhere  on  the  card, 
except  in  a body  number  where,  for  example,  1 2 
(for  body  twelve)  is  interpreted  as  1+2. 

Use  as  many  cards  of  the  above  type  as  needed  to  complete  a 
region  description,  but  leave  Columns  1-6  blank  on  all  continuation 
cards;  do  not  split  an  operator  and  its  corresponding  body  numbers 
onto  more  than  one  card. 

The  last  region  description  card  must  be  followed  by  a card 
containing  an  END  in  Columns  3,  4 and  5.  This  informs  the  code 
that  all  regions  have  been  described. 

*For  SAM-F,  the  source  region(s),  see  item  34,  must  precede 
sequentially  any  region  which  is  described  by  the  (OR)  operator. 

For  SAM-A,  the  detector  region,  (see  Section  4.4.1),  must 
precede  sequentially  any  region  which  is  described  by  the  (OR) 
operator. 


card. 

Columns 

1-2 

3-5 

6-80 
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NOTE:  If  NUSE  - 0,  Omit  Item  5 


Item  5 Character  Set  (Format  46A1 


ATABLE(N),  N«1,NUSE:  The  list  of  characters  that  are  to  be 


printed  for  each  region.  If  NUSE>47,  the  character 
is  Modulo  (46)  e.g„,  N«48  is  the  same  as  N»2.* 


The  standard  default  values  of  ATABLE  are: 


Region  Number 


1 through  9 


1 through  9 


0 (zero) 


(blank) 
Modulo  (46) 


Unless  NBRK  of  Item  4 is  >0.  Then  the  character  is  Modulo 
(NUSE-1) . 


Item  11  Option  and  Title  Card  (Format  3I3.3X 


IODT  - Cross  section  energy  banding  option 

= 0,  do  energy  banding  of  cross  sections 


= 1,  energy  banding  previously  done;  output 
must  be  available  on  logical  unit  10. 

IBEDT-  band  edit  option  (applicable  only  if  IODT-O) 

« 0,  do  not  display  energy  banded  cross  section 
arrays 

* lr  do  display  energy  banded  cross  section  arrays 

NBAND-  number  of  energy  bands  in  which  the  cross  sections 
are  to  be  processed.  Use  NBAND* 0 for  automatic 
banding  or  for  banding  previously  done  and  avail- 
able as  Tape  10.  During  tracking,  the  cross  sections 
in  only  one  band  are  in  the  computer  at  any  one  time. 
(NBAND<  49. ) Do  not  use  automatic  banding  procedure 
for  secondary  gamma  ray  problems*. 

SSDR1*  F,  no  debug  printout  from  subroutine  DR1 
= T,  debug  printout  from  subroutine  DR1 

SSDR3=  F,  no  debug  printout  from  subroutine  DR3 
= T,  debug  printout  from  subroutine  DR3 

NAME  = Title;  any  62  arbitrary  Hollerith  characters 


This  restriction  will  eventually  be  lifted 


~ ■ I 

NOTE:  If  NBAND=0,  omit  Item  12. 


Item  12  Cross  Section  Band  Limits  (Format  5E15.6) 

Enter  the  energy  limits  (in  ev)  of  each  band  starting  with 
the  highest  energy  and  proceeding  to  the  lowest  energy.  There 
must  be  NBAND+1  entries,  using  as  many  cards  as  necessary.  The 
energy  bands  must  exceed  the  energy  range  from  EHIGH  to  ECUT 
(see  Item  21  below) . 

NOTE:  If  IODT  (of  Item  11) -1, 

omit  Items  13,14  and  15. 

Item  13  Composition  Identification  (Format  2110) 

NCOMP  - total  number  of  compositions  in  the  problem. 

NNDQ  - Maximum  size  of  a band  + storage  arrays  (<NDQ) . 
Default  option:  NNDQ  = NDQ 

NOTE:  Repeat  Items  14  and  15  for 

each  (NCOMP)  composition. 

Item  14  Number  of  Elements  (Format  IIP) 

NE  - number  of  nuclides  in  present  composition  (£10). 

NOTE:  Repeat  Item  15  for  each  (NE)  nuclide 

in  present  composition. 


Item  15  Element  Card  (Format  10X,I10,E15.6) 

ID  - an  integer  which  identifies  the  nuclide  (5  decimal 
digits  ZZAAA) 

ZZ  « atomic  number 

AAA  ■ truncated  atomic  weight  for  a nuclide 

« 000  for  a naturally  occurring  mixture  of  isotopes 

CONC  - atomic  concentration  of  nuclide  in  this  composition 
in  units  of  10 2 * atoms/cm*.* 


Calculate  by  multiplying  weight  density  in  material  (grams/cm3) 

by  * - t where  A ■ non-truncated  atomic  weight  and  .6023 

■ Avogadro's  Number  in  units  of  10 2 h atoms/gm.  For  a vacuum 
description,  use  a vanishingly  small,  but  nonzero,  concentration, 
for  example,  l.E-10. 

1 


Item  21  Cutoff  Information  (Format  4E14.5) 


ECUT  Low  energy  cutoff  (ev) . Tracking  of  a particle 

is  terminated  if  its  energy  degrades  below  ECUT. 
However,  for  the  thermal  option  (see  next  entry) 
the  particle  is  not  terminated  but  is  restored, 
in  energy,  to  ETHERM..  ECUT  must  be  greater  than 
the  lowest  energy  for  which  cross  sections  are 
available. 

ETHERM  Thermal  energy  (ev) , if  a thermal  group  is 

required.  For  no  thermal  treatment,  leave  this 
entry  blank.  For  thermal  treatment,  set  ETHERM 
slightly  above  ECUT. 

IMPORTANT  NOTE: If  IODT  (item  11)  equals  0,  then 
NBAND  must  also  equal  0 on  item  11  in  order  to  use 
this  option. 

FZ  Generally,  .OKFZ^.IO.  (See  discussion 

Section  3.2.107.  At  MAGI,  FZ  » .05  is  usually 
used. 

EHIGH  High  energy  cutoff  (ev) . This  must  be  less  than 

the  highest  energy  for  which  cross  sections  are 
available. 

Item  22  Output  Energy  Bins  (Format  5E14.5)* 

These  cards  give  the  boundaries  of  the  output  energy  bins 
(ev)  used  for  the  flux  edit.  There  should  be  five  entries  per 
card  with  a total  of  (NOUT+1)  entries.  The  energies  should  be 
listed  from  high  to  low.  The  first  and  last  entries  should  be 
negative.  Any  number  of  negative  intermediate  energies  may  also 
be  used,  and  these  negative  energies  define  the  output  supergroup 
structure  (see  Section  3.2.3).  Note:  EHIGH  and  ECUT  must  be 

within  the  output  energy  bins. 


* The  complete  energy  hierarchy  for  SAM-F  is  given  in  Appendix  E. 
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Item  33  Source  Specifications  (Format  4110) 

NSR  Number  of  different  source  regions  in  the  problem. 

If  NSR-0,  an  external  source  or  a previously  gener- 
ated transmission  tape  (tape  15)  is  used  and  no  further 
source  input  is  required  after  this  card.  Skip  to  Item  40. 
If  NSR<  0 , a previously  generated  interaction  tape 
(TAPE  15)  and  a gamma  ray  production  data  tape 
(TAPE  12)  are  supplied.  Skip  to  Item  40.  For  NSR<0, 

| NSR | is  an  upper  limit  on  the  number  of  secondary 
gamma  rays  which  can  be  generated. 

IFLAG  Number  of  energies  used  in  Item  35  to  define  the 

source  spectrum.  If  IFLAG*=1 , a monoenergetic  source 
will  be  used.  If  IFLAG=0 , a built-in  Cranberg 
fission  spectrum  will  be  used  for  a neutron  problem. 
(IFLAG<150) . 

ISW  If  ISW=0,  fluxes,  will  normalize  "to  one  source  particle. 

If  ISW*1 , fluxes  will  be  normalized  to  the  total 
source  power  as  given  in  the  next  card. 

ISWW  For  external  source  (NSR=0)  only.  ISWW*0  means  dis- 

tributed source  (roust  be  isotropic).  ISWWj<0  means 
monodirectional  source  - usage  is  the  same  as  for  ISO, 
as  described  below  in  Item  34.  If  a transmission  tape 
is  used,  ISWW  is  irrelevant.  In  case  NSR<0,  ISWW  is 
automatically  set  *0. 
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NOTE:  Enter  Items  34-38  Only  if  NSR>0 


Item  34  Source  Regions  (Formats  I10,E20,,a.:riQ) 


One  card  is  required  for  each  source  region  with  a total  of 
NSR  such  cards 


ISR* 

P 


Geometrical  region  number 

Total  power  in  the  region  (source  particles  per 
unit  volume  x total  volume) . 


ISO** 


If  ISOsO,  the  source  will  be  emitted  isotropically 
(direction  may  be  biased,  however.  See  Item  29) . 


If  IS0*+1,  the  source  will  be  monodirectional , with 
the  direction  specified  on  Item  38.  No  uncollided 
estimates  are  allowable  for  point  detectors  with  ISO^O, 


If  IS0--1,  the  source  will  be  monodirectional  and  the 
code  will  not  try  to  score  uncollided  estimates  for 
small-volume  detectors.*** 


* Sources  can  be  generated  only  in  regions  consisting  of  single 
bodies.  The  body  must  be  a sphere,  right  circular  cylinder, 
box,  or  rectangular  parallelepiped.  Also,  the  source  regions 
must  be  lower  in  order  than  any  region  described  with  the  "OR" 
Combinatorial  Geometry  operator. 


**  ISO  must  be  the  same  for  all  source  regions. 


*** 


For  certain  geometries,  involving  small-volume  detectors,  no  mono- 
directional uncollided  rays  can  intercept  the  small-volume  detectors j 
Setting  IS0=-1,  will  save  computer  time,  by  instructing  the  code  not 
to  try  to  score  uncollided  fluxes.  Under  these  conditions,  IS0=+1 

is  not  incorrect,  it  is  merely  inefficient. 


o 
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Each  element  encountered  on  the  EDT  is  identified  even  if 
a particular  element  is  not  used  in  the  problem.  Thus,  if  an 
element  not  on  the  tape  is  specified  a complete  list  of  all 
elements  on  the  tape  is  printed  and  the  program  will  terminate. 

The  total  length  of  each  band  is  printed  out  and  if  re- 
quested, a printout  of  the  energy  banded  cross  section  arrays  will 
also  be  provided. 

Monte  Carlo  Data 

The  Monte  Carlo  printed  output  data  are  basically  a repeat 
of  the  input  data. 

Intermediate  Results 

All  the  parameters  of  the  first  few  source  particles  are 
printed  out.  In  number,  this  is  the  smaller  of  either  two  aggre- 
gates or  50  lines  of  printout.  The  program  will  also  print  for 
each  statistical  aggregate  the  following  items: 

Number  of  completed  histories 
Number  of  collisions 

Number  of  time  cutoffs  and  energy  degrades 
Number  of  absorptions 

Number  of  escapes  (Each  estimate  times  particle  weight) 

Total  energy  deposition  (Including  by  scattering , absorption 
and  escape  events.)  Also,  energy  of  particle  killed 
by  time  cutoff  or  energy  degrade  is  included) 

Number  of  deaths  (Including  losses  due  to  escapes,  ab- 
sorptions, time  cutoffs  and  energy  degrades.  Each 
estimate  times  particle  weight) . 

The  number  of  transmissions  and  interactions  recorded  on 
tape  14. 

Last  random  number  sequencer.  (See  Section  3.2.18) 
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preset  nor  internally  calculated,  see  input  Item  41,  then  the 
region  and  small  volume  detector  answers  must  be  divided  by  the 
region  volume.  For  point  detectors  this  is  never  necessary.) 

For  each  flux  answer  a statistical  percentage  error  is  given. 
(The  description  of  the  flux  and  statistical  percentage  error 
calculation  is  given  as  Appendix  I.)  Following  the  last  energy 
bin  the  energy-integrated  results  are  presented  for  each  response 
set.  This  is,  for  each  set,  the  sum  over  all  energy  bins  of  the 
flux  in  the  energy  bin,  times  the  energy  width  of  the  bin,  times 
the  appropriate  value  of  the  response  function  set  in  the  bin. 

3.4.3  Tape  and  Disk  File  Assignments 

The  following  are  the  file  (magnetic  tapes  or  disk)  assign- 
ments of  SAM-F.  Those  assignments  which  correspond  to  the  use  of 
various  options,  are  so  labeled.  File  numbers  refer  to  FORTRAN 
logical  numbers.  All  files  used  are  in  the  binary  mode  except  for 
file  B which  is  BCD  mode  and  contains  JBO  column  card  images. 

File  8 (Option) 

This  is  used  as  a scratch  working  file  when  automatic 
banding  is  performed.  It  is  also  used  as  the  sorted  interaction 
tape  for  secondary  gamma  ray  problems. 

File  9 

A temporary  storage  disk  file  used  by  the  BAND  subroutine. 
Subsequently,  the  file  is  also  used  for  temporary  storage  of 
particle  latents  in  the  PICK  routine,  where  it  is  referred  to,  and 
equivalenced  to,  file  18. 

File  10 


The  organized  data  file  (ODT) . The  file  contains  cross 
section  data  for  a given  problem.  The  structure  of  file  10  is 
given  in  Appendix  C.  File  10  may  be  saved  for  future  runs. 
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File  11  (Option  - may  be  omitted  if  File  10  from  a 
previous  run  is  saved.) 

The  element  data  tape  (EDT) . The  file  contains  a library 
(SAM-X  Output)  of  available  elements.  Subroutine  BAND  uses  this 
tape  to  get  the  data  for  a given  problem.  BAND  and  its  subsi- 
diary routines  are  the  only  routines  using  this  file.  The  organi- 
zation of  the  EDT  is  given  in  Appendix  A for  neutrons  and  B for 
gamma  rays. 

File  12  (Option) 

The  file  contains  secondary  gamma  ray  production  data. 

(SAM-X  Output) . 

File  13  (Option) 

This  file  contains  the  organized  secondary  gamma  ray  pro- 
duction data.  It  is  derived  from  file  12. 

File  14  (Option) 

The  interaction/transmission  file.  All  interactions  and 
transmissions  are  written  on  this  file  for  use  in  subsequent 
problems . 

File  15  (Option) 

An  external  source  file.  The  file  may  also  be  the  inter- 
action/transmission file  (File  14)  from  a previous  run,  if  it  is 

\ 

being  used  as  a supplier  of  new  source  particles. 

File  16 

The  statistical  aggregate  tape.  The  AMONTE  routine  uses 
this  file  to  record  each  completed  aggregate.  The  edit  routines 
then  process  the  file  to  obtain  the  final  flux,  and  flux-functional 
results.  Tape  must  be  saved  for  restart  option. 

File  17 

A temporary  storage  file  for  latents.  The  tape  is  used  by 
the  PICK  routine. 
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See  File  9„ 

In  addition.  Files  5 and  6 are  the  standard  input  and  output 
media,  respectively.  File  7 has  been  reserved  for  possible 
future  (or  ad  hoc)  use  as  the  standard  punch  medium. 

3.5  Variable  Core  Size  Requirement  for  the  SAM-F  Program 

At  the  present  time,  the  standard  version  of  SAM-F  loads  at 
120,000  - 130,000  octal  locations  on  a CDC  6600  machine;  the  exact 
amount  depending  on  the  given  computer  facility.  This  required 
core  allocation  can  be  reduced  or  increased  under  certain  conditions: 

(1)  The  MASTER  array  holds  the  cross  section  band  data  as 
well  as  the  output  scores  in  the  appropriate  bins.  The  size  of 
the  MASTER  array,  presently  7,000,  is  set  in  the  blank  COMMON 
statement  in  SAMF  (Main  Routine) . This  dimensioned  size  is  re- 
ferred to  as  "NDQ" • The  user  may  be  able  to  reduce  this  value 
even  further.  (Do  not  reduce  NDQ  below  4000.)  On  the  other  hand, 
the  user  may  wish  to  increase  NDQ  to  increase  the  output  scoring 
bin  capability  (see  Section  3.6)  or  to  accommodate  an  unusually 
large  bloc  of  gamma  ray  production  da<  a for  a given  nuclide. 

Note  that  the  number  of  cross  section  bands  and  thus  the 
amount  of  internal  tape  reading  and  writing  increases  with  de- 
creasing values  of  NDQ.  The  user  will  have  to  strike  his  own 
balance.  At  MAGI,  NDQ“7, 000  to  10, 000  is  usually  used. 

(2)  ISXGP  array  in  BAND,  is  presently  set  to  6630.  This  dimen- 
sioned size  is  referred  to  as  "NDD".  The  user  may  wish  to  reduce  NDD 
or  upon  bonbing  with' error  message  1205  may  have  to  increase  NDD. 

(3)  The  FPD  and  MA  arrays  are  used  in  the  geometry  routines  to 
store  body  and  region  data,  respectively.  At  the  current  dimensioned 
values  (4000  each),  about  400  bodies  and  200  regions  can  be  treated 

by  the  code.  The  user  may  wish  to  reduce  these  arrays. 
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primary  overlays „ AMONTE  can  be  reduced,  and  thus  SAM-F  reduced, 
under  the  following  conditions: 

(a)  if  no  point  detectors  are  required,  it  is  not 
necessary  to  load  ARG,  ARPREP,  FLUP  and 
TRALA. 

(b)  if  no  small-volume  detectors  are  required,  it  is 
not  necessary  to  load  FLUPV  and  TRALAV. 

(c)  if  the  thermal  diffusion  option  is  not  required, 

it  is  not  necessary  to  load  DIFFU  , DIFFUS  and  TDIFIN. 

(d)  for  primary  gamma  ray  problems  it  is  not  necessary 
to  load  ANISOT. 

3.6  Limitations  Upon  Size  of  Output  Scoring  Arrays 
The  size  of  the  MASTER  array  may  have  to  be  adjusted  (see 
Section  3.5)  to  accommodate  desired  output  scoring  array  require- 
ments „ In  addition  to  individual  array  limitations  given  in  the 
Monte  Carlo  Input  (Section  3. 4. 1.3),  the  following  limitations 
upon  NDQ  must  be  observed: 

(a)  NDQ >2  |jMAX)  (NT)  (NUMSC+NDET)]  + NUMSC 

where  MAX  = maximum  number  of  energy  bins 

in  any  supergroup 

NT  = number  of  time  bins 

NDET=  number  of  detectors 

NUMSC  = number  of  scoring  regions 

(b)  NDQ>LEGEOM  + (MAX) (NT) (NUMSC+NDET)  + 12NRMAX 

+ NRWL  + (NEWL) (NEW)  + 3NAIML  + (NUMANL) (NUMANG) 

+ 2LF0G  + 4NSR  + (JSPEC)  (NSE+1) 

+ (NUMANL+1) (NSA)  + 3 

where  LEGEOM  = storage  allocated  for  cross  sections 
(printed  out) 

NRMAX  = number  of  regions 

NRWL  =*  number  of  distinct  region  weights 

NEWL, NEW  = number  of  energy  importance  sampling 

energy  bins  and  weight  sets,  respectively 


* Sometimes  Overlay  ASORTT  is  the  largest  overlay.  ASORTT  can 
be  replaced  by  a dummy  program  if  a secondary  gamma  ray  problem 
is  not  being  run. 
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TABLE  3.3 


ERROR  STOP  MESSAGES 

STOP 


NUMBER 

ROUTINE 

OVERLAY 

MEANING 

1 

SAMF 

0,0 

Execution  went  to  completion. 

11 

GG 

0,0 

ITYPE  error.  (ITYPE<1  or 
ITYPE>10) 

21 

G1 

0,0 

Region  not  found. 

25 

SEEK 

0,0 

Vector  out  of  range. 

31 

TROPIC 

0,0 

Error  in  generating  random 
cosine  by  the  rejection  technique 

1001 

GENI 

if  o 

ITYPE  error.  (ITYPE<1  or 
ITYPE>10) 

1003 

GCHKRG 

lfO 

Too  much  data  to  continue 
geometry  checking. 

1005 

1007 

GPREE 

lfO 

Limit  of  DO  loop  exceeded. 
Limit  of  DO  loop  exceeded. 

1011 

THERM 

2,0 

Cumulative  input  data  too  great. 

1101 

AUTOBA 

2,1 

Too  many  pointers  for  element 
data. 

1103 

AUTOBA 

2f  1 

Cross  section  upper  energy 
is  below  upper  energy  of  problem. 

1105 

AUTOBA 

2f  1 

Anisotropic  continuum  inelastic 

data  given.  Coding  not  presently 
available  to  treat  this  option. 


o 
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STOP 


NUMBER 

ROUTINE 

OVERLAY 

MEANING 

1107 

AUTO  BA 

2,1 

Cross  section  lower  energy  is 
is  above  lower  energy  of  problem. 

1111 

AUTOBA 

2,1 

Elemental  data  is  missing. 

1113 

AUTOBA 

2,1 

Index  mismatch.  Memory  probably 
clobbered. 

1121 

CHITAB 

2,1 

Error  in  order  of  energies  for 
chi-tables. 

1123 

CHITAB 

2,1 

Index  mismatch.  Memory  probably 
clobbered. 

1201 

BAND 

2,2 

Too  many  composition  elements. 

1203 

BAND 

2,2 

End  of  data  tape  reached. 

1205 

BAND 

2,2 

Element  has  too  many  pointers  to 
be  processed. 

1207 

BAND 

2,2 

Error  in  processing  continuum 
inelastic  scattering. 

1211 

BAND 

2,2 

Band  of  data  is  too  large 

1221 

FILL 

2,2 

Error  in  processing  energy  tables 

1231 

LENTHS 

2,2 

Index  + 2 for  EKEBH 

1233 

LENTHS 

2,2 

Index  *1  for  E>EBLX 

1401 

INPUTD 

2,3 

Negative  region  found  for  a 
detector . 

1403 

INPUTD 

2,3 

No  room  for  input  data. 
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Either  a point  or  a volume  detector  may  be  used.  In  the  case  of  a 


volume  detector,  only  spheres,  cylinders,  rectangular  parallele- 
pipeds or  boxes  are  presently  allowed. 

(e)  Region  and  energy  dependent  importance  sampling. 

The  adjoint  transport  and  post  collision  energy  selections  are  biased 
towards  the  gamma  ray  source,  according  to  an  input  set  of  impor- 
tance values. 

(f)  Angular  Biasing  at  the  Detector. 

When  calculation  of  adjoint  currents  at  a scoring  surface  is  de- 
sired, angular  biasing  of  initial  tracks  at  the  detector  about  a 
specified  axis  is  allowed  by  means  of  specified  angular  weights. 

This  allows  aiming  toward  a surface  or  surfaces  of  interest. 

A step  by  step  description  of  the  calculation  flow  is  now  given: 

1.  Read  all  input  data  (cross  sections,  gamma  production  data, 
geometry,  importance  sampling  data,  detector  response  data). 

If  adjoint  currents  are  desired  proceed  immediately  to  Step  4. 

2.  For  a primary  gamma  ray  problem  read  in  a single  source 
particle  from  the  source  tape  and  proceed  to  Step  4.  For  a 
secondary  gamma  ray  problem  read  the  neutron  interaction  data 
for  a single  interaction.  During  the  previous  SAM-F  neutron 
problem  all  gamma  ray  producing  interactions  (absorptions  and 
inelastics)  were  written  on  an  interaction  tape.  For  each 
interaction  the  following  data  was  written  on  tape. 


(3)  neutron  collision  position 

r(3)  direction  cosines  before  scatter  (not  presently  used 

by  SAM-A  because  the  existing  coding  assumes  the 
secondary  gammas  are  generated  isotropically) 
neutron  energy  before  the  interaction 

R region  in  which  collision  occurred  (not  used  by  SAM-A) 

time  of  flight  of  the  neutron 

s anticipated  that  this  restriction  will  be  eliminated  in  the 


4. 3. 4. 2 Implementation  of  the  Biasing  Procedure 


for  Unspecified  Source  Points 


The  choice  of  the  direction  imparted  to  a ray  starting  at  the 
detector  is  done  on  the  basis  of  a biased  angular  distribution. 

A biasing  axis  is  chosen  on  input.  A set  of  polar  cosines  about 
the  biasing  axis  is  also  specified.  The  direction  of  a ray  is  then 
determined  on  the  basis  of  a biased  density  distribution  that  is 
read  in  as  input  and  the  proper  weighting  of  the  history  is  maintained 


4.3.5  Description  of  Subroutine  DOSEDT 


The  DOSEDT  routine  prints  the  detector  responses  (i.e.,  flux, 
dose,  etc.)  as  functions  of  energy  and  time  using  the  energy  mesh 
of  the  detector  response.  The  routine  is  called,  from  the  main 
program,  once  for  each  time  bin  and  also  for  the  time  integrated 
results.  For  each  time  bin,  the  collided,  uncollided  and  total 
results,  as  functions  of  energy,  are  printed  out  with  associated 
deviations.  The  energy  integrated  results  are  also  printed  out. 

Note  that  all  answers  are  normalized  to  one  primary  particle 
history  and  are  also  divided  by  the  appropriate  value  of  ^t  and  AE. 
Thus  the  results  are  on  a per  primary  source  particle  per  unit  time 
(sec)  and  per  unit  energy  (ev)  basis. 


4.3.6  Description  of  Subroutine  EDIT 


The  EDIT  subroutine  prints  the  track  length  scores  if  NG<2 


(or  the  adjoint  current  scores  for  NG*2)  and  associated  deviations 


as  functions  of  energy  and  region  azimuthal  angle  and  cosine  of  the 


polar  angle.  The  energy  mesh  used  is  EOUT  (or  ESOR)  as  described 


in  Section  4.2.5.  All  answers  are  normalized  to  one  primary  source 


particle  and  divided  by  the  appropriate  value  of  AE.  This  subroutine 


is  the  same  as  that  used  in  SAM-F  except  that  no  additional  re- 


sponse function  calculation  is  provided , since  the  detector  particle 


energies  are  selected  directly  from  the  normalized  response  function. 


4.3.7  Description  of  Subroutine  ENDET 


This  routine  is  called  from  the  main  program  to  select  the 


initial  detector  energy.  The  initial  energy  is  chosen  from  the  de- 


tector response  function  supplied  as  input. 


Item  2 - Number  of  Elements  (Format  IIP 


NE  - number  of  nuclides  in  present  composition 

Repeat  Item  3 for  each  (NE)  nuclide  in  present  composition 


Item  3 - Element  Card  (Format  10X,I10,E15.6) 


ID  - an  integer  which  identifies  the  nuclide  (5  decimal  digits 
ZZAAA) 

ZZ  = atomic  number 

AAA  = truncated  atomic  weight  for  a nuclide 

■ 000  for  a naturally  occurring  mixture  of  isotopes 

CONC  - atomic  concentration  of  nuclide  in  this  composition  in 

0 A o 

units  of  10  atoms/cnr.* 

4 - Geometry  Input 

(a)  Display  option  and  title  (Format  I10,A2,17A4) 

IPRINT  = 0,  print  out  body  and  region  data  which  follow 

= 1 , print  out  body  and  region  data  as  well  as  the  internal 
arrays  in  which  they  are  stored 
= 3,  suppress  all  geometry  printout  <• 

(b)  Body  and  Region  cards  are  the  same  as  used  by  SAM-F  and  are 


described  in  Section  3. 4. 1.1  (Items  2 and  3) 


Item  5 


N START  = Maximum  running  time  in  seconds 

**** 

NSTOP  * Total  number  of  primary  particle  histories 
NSTAT  * Number  of  primary  histories  per  statistical  group 
NRMAX  * Number  of  geometric  regions  (<_200) 

NG  * 0 , neutron  interaction  tape  supplied  for  source 
= 1,  external  gamma  ray  source  tape  supplied 
* 2,  adjoint  currents  (time- independent)  scored  in 

energy  and  angle  bins  for  tracks  leaving  scoring 
region  surfaces  abutting  the  escape  region. 

3 

* Calculate  by  multiplying  weight  density  in  material  (grams/cm  ) 
by  *6^23  where  A«non- truncated  atomic  weight  (gms)  and  .6023 
is  Avogadro's  Number  times  1024. 

**  At  present  geometry  checking  is  not  available  utilizing  SAM-A. 
Checking  of  geometry  data  should  be  done  with  SAM-F  utilizing 
the  geometry  cards  of  Section  3. 4. 1.1  (Items  1-10) 

**  Upon  completion  of  each  statistical  aggregate  the  code  estimates 
the  total  (real)  running  time  at  the  end  of  two  additional  aggre 
gates.  If  this  time  exceeds  NSTART  (seconds)  no  new  aggregates 
will  be  started  and  control  is  switched  to  the  edit  routines. 

**  Refers  to  primary  neutron  histories  for  secondary  gamma  ray 
problems. 


MT  = Number  of  time  bins  (<_25 );  for  no  time  dependence  use  MT= 
(No  time  dependence  is  currently  allowed  for  NG»2  option) 

NOUT  « Number  of  output  energy  bins  (<_25) 

MRS  * Number  of  scoring  regions 


NRWL  = Number  of  distinct  regions  weights  (<  101) 


IESC  = Escape  region  number 


IBG  = 0,  no  debug  printout 


If  NG*=0,  MUL  is  the  number  of  times  each  neutron  interaction 
is  to  be  used.  (If  a zero  is  entered,  the  code  assumes 
MUL=1) 


Item  6 - Detector  Parameters  (Format  3B14.6,I4) 


XD  * Detector  coordinates  (x,y,z);  (for  point  detector  only) 

* 

IRDET=  Detector  region  number  (use  negative  number  for  point 
detector,  or  a positive  number  for  volume  detector) 


Item  7 - Output  Gamma  Energy  Mesh  (Format  5E14.5) 


The  mesh  at  which  the  track  length  scores  will  oe  tallied.  Enter 


NOUT+1  numbers  from  high  to  low  energy.  Use  as  many  cards  as  required 


Item  8 - Output  Time  Bin  Mesh  (Format  5E14.5) 


Time  mesh  (seconds)  at  which  detector  scores  will  be  tallied 


Enter  MT+1  numbers  from  longest  to  shortest  times.  Omit  this  set  if 


Item  9 - Region  Weights  (Format  514.6) 


Enter  the  NRWL  region  weights  to  be  used  in  the  problem.  The 
weights  need  not  be  entered  monotonically  by  value  but  their  order 
determines  the  region  weight  numbers  (i.e.,  entry  one  is  weight  #1 
etc.) . 


The  volume  detector  must  be  a region  restricted  to  a single  body. 
Either  a sphere,  right  circular  cylinder,  box,  or  rectangular 
parallelepiped  is  allowed.  Furthermore,  no  regions  employing  the 
(OR)  operator,  see  Section  3. 2. 1.1,  may  sequentially  precede  this 
detector  region. 


o 


NOTE:  Enter  Items  17-24  only 

if  NG-2  (See  Item  5) 

Item  17  - Scoring  Parameters  (Format  3110) 

MWSR  - Number  of  polar  angle  cosine  bins  <_34 
MSOR  - Number  of  scoring  energy  bins  <_34 
MASR  - Number  of  equally  spaced  azimuthal  bins  <_34 

Item  18  - Scoring  Polar  Angle  Cosines  (5E14.5) 

Polar  cosine  mesh.  Enter  MWSR+1  cosines,  in  descending  order, 
from  +1.  to  -1. 

Item  19  - Scoring  Energies  (5E14.5) 

Scoring  energy  mesh  (ev) . Enter  MSOR+1  numbers,  in  descending 
order . 

Item  20  - Scoring  Axis  Direction  Cosines  (3E14.5) 

Specification  of  direction  cosines  of  the  axis  for  polar  cosine 
mesh  as  specified  in  Item  18.  Enter  x,y,z  direction  cosines 
in  order. 

Item  21  - Biasing  Axis  Direction  Cosines  (3E14.5) 

Specification  of  aiming  axis  about  which  is  specified  a polar 
cosine  mesh  for  biased  selection  of  the  direction  out  of  the 
detector  (Item  23).  Enter  x,y,z  direction  cosines  in  order. 

Item  22  - Number  of  Biasing  Polar  Cosines  Bins  (IIP) 

Enter  number,  NB,  of  such  bins.  <20 

Item  23  - Biasing  Polar  Cosine  Mesh  Points  (5E14.6) 

Enter  NB+1  cosines,  in  descending  order,  from  +1.  to  -1. 

Item  24  - Biasing  Polar  Cosine  Direction  Weights  (5E14.6) 

Enter  NB  weights. 

4.4.2  Description  of  Output 

The  SAM-A  output  falls  into  the  following  general  categories 
(in  the  order  in  which  they  are  printed) . 

a.  Cross  section  and  composition  information, 

b.  Geometry  information, 

c.  Problem  input  data, 

d.  End  of  aggregate  tally, 

e.  Track  length  or  adjoint  current  scores, 

f.  Gamma  ray  collision  tally, 

g.  Detector  edit.  (NQ< 2) 
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After  the  final  aggregate  the  total  Monte  Carlo  running  time  is 
printed. 
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e.  Track  Length  or  Adjoint  Current  Scores 

NQc 2 - For  each  region  designated  as  a scoring  region,  (input 
parameter) , the  track  length  scores  as  function  of  the  track  length 
energy  mesh  (also  input)  will  be  printed.  The  displayed  results  have 
been  divided  by  AE  and  by  the  number  of  primary  histories.  This  edit 
is  time-integrated  and  is,  in  general,  the  same  as  used  for  track  length 
scores  in  the  SAM-F  program.  Note  that  the  answers  are  not  divided 
by  the  region  volumes. 

In  addition  to  the  energy  dependent  scores,  the  energy  integrated 
scores  are  also  given  for  each  region. 

NG=2 

For  a single  azimuthal  bin  and  a single  scoring  region  the  adjoint 
current  computed  from  the  track  length  scores  are  given  as  functions 
of  the  track  length  energy  mesh  and  the  cosine  of  the  polar  angle. 

Both  energy  dependent  and  energy  integrated  scores  are  given.  The  ad- 
joint current  printout  is  given  for  each  scoring  region.  When  more 
than  one  azimuthal  bin  is  given  the  adjoint  current  is  printed  out 
only  for  the  first  bin  in  order  to  shorten  the  otherwise  voluminous 
output . 


f.  Gamma  Ray  Collision  Edit 

■ 

The  number  of  gamma  ray  collisions  in  each  geometric  region,  (as 
a function  of  energy,  using  the  track  length  energy  mesh) , is  printed. 
Region  numbers  are  given  across  the  page  and  energy  bin  bounds  are 
given  down  the  page. 

g.  Detector  Response  Edit  (NQ< 2) 

The  detector  response  edit  give  the  uncollided,  collided  and  total 
detector  responses  as  functions  of  energy  and  time.  The  results  are 
scored  in  bins  using  the  energy  mesh  for  the  response  function  data 
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and  the  input-specified  time  mesh.  The  units  of  the  response  corres- 
pond to  the  detector  response  function  supplied  as  input.  Hence,  if 
the  input  response  function  was  flux-to-dose  conversion  factors,  the 
final  results  are  doses.  If  a unity  response  function  was  supplied,  ! 
the  results  are  number  fluxes.  If  the  average  energy  of  a scoring 
energy  bin  was  entered,  the  results  are  energy  fluxes,  etc. 

Each  time  bin  is  edited  separately  and  the  energy  integrated  re- 
sponse is  given  for  each  time  bin.  In  addition,  a time-integrated  re- 
sponse is  given  after  the  time-dependent  results  for  each  time  bin 
have  been  printed. 

4.4.3  Tape  and  Disk  File  Assignments 

The  following  are  the  file  (magnetic  tapes  or  disk)  assignments 
of  SAM-A.  File  numbers  refer  to  FORTRAN  logical  numbers. 

File  7 (Option) 

The  gamma  ray  production  cross  section  file. 

File  11 

The  gamma  ray  cross  section  tape. 

File  14 

If  NG<  2 it  is  the  neutron  interaction  tape  or  the  primary  gamma 
source  tape.  When  NG=2  the  adjoint  currents  are  written  on  tape  in 
the  form  of  BCD  records. 

In  addition.  Files  5 and  6 are  the  standard  input  and  output 
media,  respectively. 

When  NG-2  the  interactions  or  gamma  ray  sources  are  written  on 
the  tape  in  blocks  of  35  particles  per  record  with  the  following  write 
statement. 

WRITE (14)  ( (B(I,J) ,1=1,14) ,J*1,35)  | 

Thus,  each  column  of  ”B”  is  a 14  word  particle  record. 

The  format  of  each  14  word  particle  record  on  File  14  is  shown 
below. 

a 

- 
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Rearranging  equation  (G.2)  and  squaring  both  sides 


Substituting  equations  (G.la,b,c)  into  equation  (G.3)  yields 
after  canceling  common  terms  and  rearranging 


where  the  coefficients  of  the  quadratic  equation  are  given  by 


4 (1-u  ) Ef + (A+l)  r +4 (1-w  +A) rE 


4 (u+A-DrE,  +2  (2u>  +A  -1)  r +4ru>  (E.  +r)  /w  -1+A 


2 (A+l)  r +8 (l-w  )ET+8 (1-w  +A)rE 


where  we  have  retained  the  root  which  yields  T»1  for  u»l  (recall 
that  the  additional  root  was  introduced  by  squaring  equation  (G.2)) 


k* 


tVuX«i.5 

|x-xkl2  lx-x£l2 


) - 


Vi 

| X"X.  | | X-XjJ 


the  COT  is  satisfied,  and  an  MDC  exists. 

Although  references  to  "sphere  overlaps"  (see  algorithm  (0-2)) 
and  "multiple  intersections"  (see  algorithm  (0-3))  also  imply  the 
existence  of  an  "angle"  MDC,  these  tests  comprise  a sub-set  of  the 
COT.  Finally,  only  a sphere  overlap  is  a significant  MDC  for 
position  reselection. 

4.3  Resolution  of  an  MDC 

An  MDC  may  be  resolved  by  Russian  roulette,  in  which  one 
("live")  detector  is  chosen  for  sampling  purposes.  The  implication 
of  this  selection  depends  on  the  particular  stage  of  the  particle 
history  at  which  it  occurs. 

4.3.1  Position  Checking  (Source) 

After  the  selection  of  a source  position  at  the  start  of  a 
history,  the  resolution  of  an  MDC  is  recorded  by  the  appropriate 
setting  of  the  live  detector  "position  flag",  KDLIV.  The  relevant 
algorithms  are  best  summarized  by  a logical  flow  chart  (Figure  0-2). 
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4.3.4  Interpretation  of  KDLIV 

As  a result  of  positon  checking,  both  for  a source  and  a 
collision,  the  setting  of  KDLIV  has  the  following  implications: 

(a)  The  case  KDLIVM)  implies  that  there  was  no  MDC. 

(b)  For  KDLIV<0,  detector  number  ( -KDLIV)  is  used  for 
sampling  purposes,  and  a special  weight  adjustment 
allows  bounded  estimation  for  all  detectors.  This 
weight  adjustment  is  given  by 

n b.  /I  b 
k m=l  m 

where 

k - | KDLIV | 

V * 1 ' Rm2/'  I x_xm  1 2 * 

(c)  The  value  of  KDLIV  influences  the  direction  checking 
preceding  angle  reselection,  as  described  in  Section  4.3.2  above. 
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The  data  sheets  which  follow,  provide  the  complete  input 
description  for  Combinatorial  Geometry  Example  #5,  as  shown  in 
Section  3. 2.1„3. 


SPECIAL  INSTRUCTIONS  IDENTIFICATION 


SECTION  6;  INDEX  (BY  PROGRAM) 


Collision  Mechanics  187-192 
Elastic  Scattering  187-188 
Inelastic  Scattering  188-189 
Continuum  189-190 
NCDB  Parameter  192 


BCDEAN  (See  Output  Format) 

Core  Size  Requirements  46-48,  112-113 
Cross  Sections 


ENDF  Data  Files  15-20,  33,  et  passim 
GAMMA  (See  Gamma  Ray  Cross  Section  Data) 
Gamma  Ray  Cross  Section  Data  17,45,95 
Subroutine  Descriptions  95 


COMMON  Arrays  212-216 

Core  Size  Requirements  201,  246-248 

Cross  Sections 

Gamma  Ray  Element  Data  Tape  319-320 
Gamma  Ray  Production  Data  Tape  323-325 
Neutron  Element  Data  Tape  305-317 
Organized  Data  Tape  321-322 
Processed  Data  Tapes  Currently 
Available  96b 
Thermal  186-187 
Use  in  Monte  Carlo  191-192 


Gamma  Ray  Production  Data  16-17 
79-94 

Subroutine  Descriptions  79-94 


General  Descriptions  15-17,  46-48 
Input  Descriptions 
BCDEAN  103-108 

General  Discussion  97,  109-110 
Main  (SAMX)  97-98 
NUTRON  98-103 

PEND  103-104 

WEED  105-106 


Energy  Hierarchy  327 
Error  STOP  Messages  248-252 
Escape  Regions  162 
Flux-At-A-Point  (See  Scoring) 

Flux- In- A- Small  Volume  (See  Scoring) 
Fourteen  Word  Particle  Array  207 
Geometry  (See  Combinatorial  Geometry) 
General  Descriptions  115-117,  163-168 
IDET  Parameter  203,207 
Importance  Sampling  150-155,  182-184, 
205,  329-338,  397-400 
Input  Descriptions  217-241 
Geometry  218-227 
Cross  Sections  228-229 
Monte  Carlo  230-241 


NUTRON  (£  o Neutron  Cross  Section  Data) 
Output,  and  Output  Formats  17,  45,  78, 

96,  305-317,  319-320,  323-325, 

.t  71-373 

lay  Structure  47 

Pi  (see  Gamma  Ray  Production  Data) 

Pr<  ^ussed  Cross  Section  Data  96a-96b 
in-,  .ction  Types  33 

Resonance  Region  (See  Cross  Sections) 
Secondary  Angular  Distributions  36-39, 
64-74 

secondary  Energy  Distributions  39-42,74-78 
Tape  and  Disk  File  Assignments  110-112 
WEED  (See  Gamma  Ray  Production  Data) 


Interaction  File  160-161,  179-180,  206 

Interactions  182-185,  191-192 

J12345  Parameters  203,207 

Latents  154-155,  176-177,  185,  201-203 

MASTER/ASTER  Array  343-345 

NCDB  Parameter  192 

NDD  Parameter  246 

NDQ  Parameter  201,  246-248 

Output  Descriptions  241-244 

Aggregate  (Intermediate)  Tally  209-210, 
242 


i . lOGRAM  SAM-F 


ending.  Supergroups,  Superbins  142-143 
176-177,  180-181,  201-203 
Collision  Events  (See  Interactions) 


OVERLAY  Structure  163-168,  211 
Random  Number  Sequencer  162,  230 
Relativistic  Effects  190,  339-342 


PROGRAM  SAM-F  (Continued) 


Response  Functions  160,  241 
Restart  Option  162,  230 
Scoring  116,  148-149,  347 

Flux-At-A-Point  178-179,  193-195, 

210,  375-391 

Flux-In-A-Small  Volume  196-200, 

211,  236 

Track  Length  116,  347 

Secondary  Gamma  Ray  Production 
179-180,  206 

Source  Capabilities  144-146,  179-180, 
204-209,  236 

Statistics  158,  230,  347-348 
Subroutine  Descriptions  (Detailed) 

176-210 

AMONTE  176-177?  ARG  178-179? 

ASORTT  179-180?  BAND  180-181? 

CARLO  181-185?  CARSCA  185-186? 

DXFIN  186-187?  DR3  187-192? 

FLUP  193-195?  FLUPV  196-200? 

MAIN  (SAMF)  201?  PICK  201-203? 

SEEK  204?  SOUCAL  204-209? 

SOUGAM  206?  SOUGEN  206-209? 

SOUPIC  206-209?  TALLY  209-210? 

TRALA  210?  TRALAV  210 

Subroutine  Descriptions  (General) 

169-175 

Superbins  (See  Banding,  etc.) 

Supergroups  (See  Banding,  etc.) 

Tally  (by  aggregate)  242-243 

Tape  and  Disk  File  Assignments  244-246 

Thermal  Neutrons  156-157 

Diffusion  Option  157,  186-187,  239 
Tracking  181-185 

Transmission  Regions  161,  401-404 
Volume  Computation  159 
WC  Parameter  207 

PROGRAM  SAM-A 

Collisions  (See  Tracking  and  Collisions) 
Combinatorial  Geometry  (See  SAM-F  Section) 
Cross  Sections  258,  265,  278-279 
Gamma  Ray  Element  Data  Tape  319-320 
Gamma  Ray  Production  Data  Tape  323-325 
Organized  Data  Tape  321-322 

Detector  Particle  Specification  261,  270, 
279-283,  290-292,  295,  299 


Energy  Meshes  and  Hierarchy  259-260,  262 
Escape  Region  264 

Fourteen  Word  Particle  Array  255-256, 
289-290,  301 

General  Descriptions  253-257,  264-271 
Geometry  (See  SAM-F  Section, 

Combinatorial  Geoeietry) 

Importance  Sampling  263,  279-282,  349-362, 
Also  See  Sam-F  Section,  Importance 
Sampling 

Input  Descriptions  293-297 
Cross  Sections  293-294 
Geometry  294,  218-227 
Monte  Carlo  294-297 

Interaction  File  255-256,  289-290 
Output  Descriptions  297-302 

Aggregate  (Intermediate) Tally  298 

Scoring  284-289,  292-293,  299,  405-413 
Source  Capabilities  261,  289-290 
Statistics  263,  294,  347-348 
Subroutine  Descriptions  278-293 
DATORG  278-279;  DIREC  279; 

DIRSEL  279-282;  DGSEDT  282; 

EDIT  283;  ENDET  283;  ESTMT  284-289; 
GASSUP  289-290;  POSGEN  290-292? 

SET  292?  SEEK  292?  UNCOL  292-293 

Tape  and  Disk  File  Assignments  300 
Time  Dependence  262 

Tracking  and  Collisions  270-277,  349-362, 
397-400 
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